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Abstract 

We show that in weakly non-ideal plasmas, like the solar interior, both non- 
extensivity and quantum uncertainty (a la Galitskii and Yakimets) should 
be taken into account to derive equilibrium ion distribution functions and to 
estimate nuclear reaction rates and solar neutrino fluxes. 
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It is well known that Boltzmann-Gibbs statistics does not describe correctly the statistical 
behavior of particles with long-range interaction, non-Markovian memory and multifractal 
boundary conditions [Q,^ . The treatment in the frame of the extensive description is based 
on short-range interaction, on neglecting surface effects and on fluctuations of thermody- 
namic quantities relatively too small to be observed. 

The solar core plasma represents a system of particles that, for many reasons, has an 
equilibrium distribution deviating very slightly from the Maxwell-Boltzmann (MB) distri- 
bution. In fact, the mean Coulomb energy is not much smaller than the thermal kinetic 
energy; the Debye length is of the order of the inter-ionic distance a^; it is not possible to 
clearly separate collective and individual degrees of freedom; the presence of the scales of 
energies of the same size produces deviations from the standard statistics [Q. 

In the region of plasma parameter F = Z'^e^/{aikT) k, 0.1, to which the solar core 
plasma belongs, the Debye character of the screening is maintained but strong interactions at 
distances of the order of n~^^^ [n is the particle density) become the main interactions while 
involving multiple particles. The reaction time to build up screening after a hard collision is 
the inverse of the plasma frequency; the collision time is comparable to it. Therefore many 
collisions are necessary before particles loose memory of the initial state and the scattering 
process can not be considered Markovian, non-extensive effects should be included 0. Tsallis 
non-extensive thermostatistics is actually a well-known tool successfully applied to many 
different physical problems and can account quite naturally of all the above mentioned 
features of the solar core plasma by means of the entropic parameter q 

Recently, in Ref.s |^,^ the Green's function technique has been used to examine tunneling 
reaction probability and the role of the many-body effects (or quantum uncertainty effect, 
following Galitskii and Yakimets 0) on the rates of nuclear reactions in the laboratory 
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fusion plasma and in the weakly non-ideal sun core. The implications of such effects in the 
evaluation of the solar neutrino fluxes are discussed in Ref. P] (for an updated reading on 
the solar neutrino problem one can refer to Ref.s [p|-p!0[] and to the recent Ref. |]1T| where 



neutrino decay solution is reexamined in the light of the Super Kamiokande data). 

The aim of this paper is twofold: 1) we find a microscopic interpretation of the non- 
extensive parameter q in the solar core (in other words, the entity of deviation of the 
Maxwellian distribution of the ionic equilibrium distribution we use in this work) in terms of 
the plasma parameter and the ion correlation parameter; 2) we study both the non-extensive 
(g) and the quantum uncertainty {Q) effects [^[Tj, in the calculation of the solar reaction 
rates and in the solar neutrino fluxes. The two effects should be simultaneously active if the 
collision frequency is comparable with the plasma frequency; this requirement is fulfilled in 
the solar core. 

Electrical microfields and validity of the Tsallis distribution in the solar core - Basic 
quantities influencing the microscopic processes and the nuclear rates are the electrical mi- 
crofields acting on the nuclei in the plasma: the time-spatial fluctuations in the particles 
positions produce specific fluctuations of the microscopic electric field in a given point of 
the plasma. The rates are changed respect to the standard ones because of the presence of 
microfield distributions (with energy density of the order of 10~^^ MeV/fm^) which modi- 
fies, under particular conditions, the ion distribution tail. The influence of microfields on 
the rates of nuclear tunneling reactions has been widely studied since the early work of 



Holtsmark and later, for instance, in the Ref.s [|T^-|T6|. 

In the solar core, the effects of random electric microfields are of crucial importance. 
These microfields have in general long-time and long-range correlations, can generate anoma- 
lous diffusion and may be decomposed in three main components: 1) A slow varying com- 
ponent due to plasma oscillations. The particles see this component as an almost constant 
external mean field over several collisions. 2) A fast random component related to the diffu- 
sive cross section {ad ~ 1/^); this component does not affect the distribution that remains 
Maxwellian. 3) The third component is related to a short-range two-body strong Coulomb 
effective interaction. As we will see later, this is the component that alters the distribution. 

Let us derive the ion distribution function to be used to calculate tunneling reactions 
rates in solar plasma, when an electric microfield distribution is present. The equilibrium 
distribution we are deriving differs from the MB distribution, if particular conditions are 
fulfilled. In fact, the presence of the electric microfield average energy density, {E"^), modifies 
the stationary solution of the Fokker-Planck equation and the ion equilibrium distribution 



can be written as [|T7| 



3 xmi/^ 



where u is the total collision frequency, x the elastic energy-transfer coefficient between two 
particles of the plasma, x = 2mim2/(mi + m2)^, m the reduced mass, T the temperature, 
C the normalization constant. 



Defining a critical field Ec = z/y 3xm/i;T/2e^, we can see from Eq.(|l]) that ii E E^ the 
distribution is Maxwellian whatever be the value of the frequency u, if E ^ E^ and u is 
not a constant but depends on v, the distribution is a Druyvesteyn like distribution |T^. In 
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the solar core being E not too larger than Ec, the distribution of Eq.(|l|) differs slightly from 
the Maxwellian but such small deviation is quite important in the evaluation of the nuclear 
rates. 

The condition that z/ be a function of velocity v is verified by the fact that the elastic collision 
cross section is a = + ctq, where oc 1/f is the elastic diffusion cross section and o"o is 
the enforced Coulomb cross section. The total frequency u = (avn) satisfies the relation 
i/^ = i^J + z/q without interference terms. 



The expression of ctq is due to Ichimaru [T^ which developed a strict enforcement of the 
Wigner-Seitz ion sphere model yielding the elastic cross section ctq = 27r(aa)^ where a is the 
inter-particle distance, a is a one dimensional parameter related to the probability that the 
nearest neighbor be at a distance R and therefore related to the pair-correlation function 
giR,t). 

The explicit expression of the equilibrium distribution ([^) for the solar interior can be 
written as a function of the kinetic energy 



/(cp) = A^exp 



where N is the normalization constant and 



9 n^{kTf 



(2) 



(3) 



0"n 




In solar interior ip ^ 1 {if ^ 1), therefore the equilibrium distribution function contain- 
ing random microfields with collision frequency depending on the velocity, is given by the 
Maxwellian distribution times the factor exp[—6{ep/kT)'^]. It represents, if we recall 6 as 
5 = (1 — g)/2, the approximation of the Tsallis distribution when q ^ 1 [0. It is remarkable 
that such distribution has been postulated ad hoc by Clayton more than twenty years ago 



in the solution of the solar neutrino problem |19 



By means of Eq.s(|^)-(^) we have established a strict connection between the entity of the 
deviation from the MB in the solar interior, the 6 (q) parameter, and the elastic diffusion and 
Coulomb cross sections. Expliciting such cross sections in terms of the Ichimaru parameter 
a and the plasma parameter F, after straightforward calculations, we obtain that, in the 
small correction limit relevant to the solar core (v? ~ 1), the 6 parameter can be written as 

|5|^-f^ = 12a^r2«l. (5) 

Eq.(^ is a crucial result of this paper because it establishes that the presence of electric 
microfields implies a deviation from the MB distribution, the entity depending on the value 
of the plasma parameter and the collision frequency by means of the a parameter. As we 
will see, very small deviations from the Maxwellian tail can produce strong deviations of the 
reaction rates from their standard values. 
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Let us remark that we do not discuss the electron screening effect pO| . We can reahze that 
at the temperature, elemental densities and plasma parameter of the sun core, the rates of the 
reactions could be modified of only few per cent and even less by the screening factor /q. Its 
value can be taken equal to unity because its effect on the rates is negligible if compared to the 
depletions (or enhancements) of the rates due to other effects (electromagnetic fluctuations 
in a plasma) responsible of the equilibrium distribution function we are considering and 
evaluating in this work. 

The quantum uncertainty effect - Galitskii and Yakimets showed that in dense or low 
temperature plasma, due to quantum uncertainty effect, the particle distribution function 



over momenta acquires a power-like tail even under conditions of thermodynamic equilibrium 

p 



f{e,p)de = fM{p) + ^e^/'\ (6) 



where /Af(p) is the equilibrium MB distribution over momenta. 

The distribution function of Eq.(|^) is an approximation of the correct solution, not valid 
for hv/kT ^ 1 (non-ideal plasma of high density and low temperature), and is correct in 
the case of solar core hu/kT ~ 0.1 (we remember, as already mentioned, that u ~ ujpi). 
At sufficiently high values of p, the magnitude of the tail can be many times greater than 
the Maxwellian distribution. Since is the momentum rather than the energy that enters 
into the scattering amplitude in the gas approximation, fusion reaction rates in equilibrium 
systems can be calculated by averaging the reaction cross section cr(ep) over the momentum 
distribution function, rather than the energy distribution 0. 

Let us briefly summarize the results obtained in Ref. 0] taking into account the uncer- 
tainty quantum effects Q alone. We indicate the corrections derived respect to the standard 
results based on Maxwellian distributions. 

Although the pp reaction rate should be corrected by a factor ru ~ 3.5 ■ 10"'^, (at 
r = 0.0759 i?©) nevertheless the solar luminosity is maintained unchanged respect to the 
observed value because: a) the ^He-^He rate should be increased by a factor r^^ ~ 4.5 ■ 10^, 
b) the ^He-^He rate should be increased by a factor r^^ ^ 3.1 • 10^, at r = 0.00759 i?© and 

~ 5.9 • 10^, at r ~ 0, c) the ^He density is unchanged, d) ^He density is decreased by 10^, 
e) electronic capture (that is not a tunneling reaction) rate by ^Be is unchanged, f) p-^Be 
rate is increased by a factor ~ 4 ■ 10^, at r ~ 0, g) the density of is increased by 8 ■ 10^, 
h) other reaction rates to burn ®B are increased by ~ 10'^°, i) the flux of ''Be neutrinos is 
decreased by a factor ^ 50. 

The authors argue that because the mean average molecular weight A is left unchanged 
then, the sound speed does not change respect to the standard prediction. This should be 
proved, because locally in interior shells, where several reactions are mostly active, different 
density (from standard) might induce quite different speed of sound and helioseismology 
constraints can be not fulfilled. The aim of Ref. was to check if the calculated quantum 
uncertainty Q rates contradict solar data and found it do not. Solar interior conditions 
require the inclusion of the Q corrections. Nevertheless, because results are quite different 
from solar standard model (SSM) results we will see that only within non-extensive ther- 
mostatistics these corrections can be adequately arranged to agree with the experimental 
results. 
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Nuclear reaction rates with non-extensive and quantum uncertainty corrections - We 
know that a solution which coincides with a non-extensive distribution can be obtained only 
if the electric microfield E is greater than a certain critical value i?c, and if the frequency 
of collision depends on v. This condition is verified if the Coulomb collisional cross section, 
which describes the strong part of the collision, is a constant: ctq = 27r(aa)^ and is due to 
the random contribution of the electric microfield. The above assumptions do not modify 
the treatment of Galitskii and Yakimets (they consider the Coulomb interaction, but the 
random part of it is not taken into account), rather are responsible of some changements in 
the final rates. 

The non-extensive distribution function, given by Eq.(|^), modifies the reaction rate Rq 
and can be written as 

i?, = i?Me-'^- , (7) 

where Rm are the Maxwellian rates, 5 = (1 — q)/'2., ■jij = (eg/fcT)^, eg = 
5.64:{ZiZj AiAj / (Ai + Aj)Tc/Ty^^, is the most effective energy and Tc = 1.36 keV is the 
temperature at the center of the Sun. 

If we consider only the quantum effect Q, the correction rij to the Maxwellian rate 
{Rq = RuTij) can be written as [|| 



87r3/2 ^hT\ m, / r,^- 

rricoii = ■ , rur = ■ , (9) 

+ mb mi + rnj 



where ruh is the mass of the background particles colliding with tunneling particles. 

Including the Tsallis non-extensive effects, the total rate (i.e. including both the q and 
Q effects) is 

R = Rq + Rqe-'-'h = (e-^^- + r,^e-'^h) , (11) 

where 7* = {tq/kTY and eg is the effective energy of the quantum correction. Consider- 
ing the appropriate expression of the collision frequency v in terms of ctq and 0"^, we can 
verify that approximately holds the numerical relation: eg ~ 3eo. Although the values r^j 
calculated in Ref. [^] are of order 10® -i- 10^, the factor e~^'^* has a large exponent and can 
strongly suppress the enhancement given by r^j. In fact, the quantum uncertainty correc- 
tions calculated on the basis of Tsallis distribution rije~^^* are of the same order of e'^"^ and 
even less (if 5 ~ 10^^ ^ 10^^, as we can deduce from Eq.(^). Solar conditions admit small 
non-extensive statistics effects in the equilibrium distribution, showing very small deviations 
in the Maxwellian tail. By using a numerical code based on a complete evolutionary stellar 
model, it has been verified that the consequence of this statement are compatible to the 
experimental results on neutrino fluxes |2T| . 
Let us now define the following quantities: 
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and /ce? = -2^17, ke^7 = ykn, where kij and kfj are respectively the modified and MaxweUian 
rates; the value of the constant z = 227 can be found, for instance, in Ref. ||22||. 



We can take the usual time-dependent equations of the rates with the solar luminosity 
constraint (as reported for instance in Ref. 0) to derive the steady state solutions for the 
elemental densities. Using the expressions of the rates of Eg. (pi]) with quantum uncertainty 
and non-extensive corrections, the following set of equations can be derived: 

A=- = e-'^'^'' < nye-'i^^i^ , --7 = 1, 

B rn z A 



n^g7 n. 



1 + (e-'5i77ir + e-'^i^Ti*7)/(2z)] ■ 



A reasonable evaluation of a gives: a = 0.55, with F ~ 0.1 we obtain q = 0.990 {6 = 0.005) 
for all components (see the Conclusions for a comment on the uncertainty of these values). 

If we assume n^/n^^ 3 • 10^'^ (let us remark that density is not changed alone, but all 
the tunneling rates are changed consistently solving the set of the mentioned equations) we 
obtain 

$(5e^) 1 $(5) 1 



$^(5e7) 50 ' $^^(5) 2 ' 
Gallium = 81 SNU , Chlorine = 2.8 SNU , 

$(pp) and luminosity are practically unchanged respect to the SSM values. 

The CNO rates that are strongly enhanced by the quantum uncertainty effect are re- 
markably reduced by the factor e"'''^*. We have introduced the assumption that is of 
the order of about 10~^ reduced respect SSM and consistent with the set of time-dependent 
equations of the rates. This value is within the constraints actually imposed by helioseismol- 
ogy because in the region r/Rgun < 0.2 the value of can be submitted to a large variability 
|23| . The results presented, that are in very good agreement with the experimental data 
, do not pretend to demonstrate that the solar neutrino problem is solved, but rather 



that a path to solution along this procedure can be found. 

Conclusions - In a weakly non-ideal plasma, like the solar interior, ions and electrons 
behave like quasi-particles and, due to the quantum uncertainty, momentum and energy of 
colliding particles are independent variables. This many-body effect is responsible of an ad- 
ditional power-law tail to the Maxwellian momentum equilibrium distribution, as discussed 
by Galitskii and Yakimets 0, while the energy particle distribution remains Maxwellian. 
The central and more important quantity that enters into the quantum correction of the 
distribution is the elastic collision frequency which, in these weakly non-ideal plasmas, has 
a value of the same order as the collective plasma frequency. Being such a value of colli- 
sion frequency in systems of particles with long-range interactions related also to memory 
effects with long-time tails, the equilibrium non-extensive statistics appears to be the most 
suitable statistics to use and it can take into account, in a very natural way, both quantum 
and non-extensive effects. 
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Savchenko |5| has shown that the rigorous procedure of averaging of tunnehng probabihty, 
based on Green function technique, gives results nearly coincident with the results obtained 
by the simple averaging. Therefore, we can show that, by using the modified distribution 
function, Eq.(0), (which corresponds to the equilibrium Tsallis distribution when deviations 
from standard statistics are small) instead of the MB one, the strong increase of the nuclear 
reaction rates caused by quantum uncertainty effects alone is greatly (or fully) reduced. We 
can confirm the validity of the Tsallis statistics to describe weakly non-ideal plasmas such 
as the solar core and of its use to calculate solar neutrino fluxes towards a closer agreement 
to the experimental measurements than standard calculations. In fact, the calculated fluxes 
are in good agreement with the experimental results using the reaction rates of Eq. (|1T]) with 
quantum uncertainty and non-extensive corrections and a ^He density (in the core) 10~^ 
times the value usually taken in SSM. This figure is consistent with helioseismology con- 
straints. Finally, we have justified the entropic parameter q (5) of the non-extensive statistics 
because we have related this quantity to the electric microfields distribution in the plasma 
(or explicitly, to the plasma parameter F and to a, the ion correlation function). Therefore 
the quantity q cannot be considered a free parameter. However, the evaluation of q (or 5) by 
means of Eq.(^) is not absent of some uncertainty because the quantity a (5 depends on the 
fourth power of a) is still to be deeply analyzed and evaluated, in spite of the many efforts 
spent in the past on space-time correlation functions of ions in stellar plasma. Our results 
are, of course, not conclusive on the problem of solar neutrinos, rather they represent an 
indication that, possibly, an astrophysical solution can be found once nuclear reactions in 
weakly non-ideal astrophysical plasma be fully understood or measured in laboratories with 
ad hoc experimental facilities. 

This work is supported in part by MURST. 
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